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Abstract 

The electronic excitation population and coherence dynamics in the chromophores of the 
photosynthetic light harvesting complex 2 (LH2) B850 ring from purple bacteria {Rhodopseu- 
domonas acidophila) have been studied theoretically at both physiological and cryogenic 
temperatures. Similar to the well-studied Fenna-Matthews-Olson (FMO) protein, oscilla- 
tions of the excitation population and coherence in the site basis are observed in LH2 by 
using a scaled hierarchical equation of motion (HEOM) approach. However, this oscillation 
time (300 fs) is much shorter compared to the FMO protein (650 fs) at cryogenic tempera- 
ture. Both environment and high temperature are found to enhance the propagation speed 
of the exciton wave packet yet they shorten the coherence time and suppress the oscilla- 
tion amplitude of coherence and the population. Our calculations show that a long-lived 
coherence between chromophore electronic excited states can exist in such a noisy biological 
environment. 
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I. INTRODUCTION 



Photosynthesis is the single most crucial biochemical process in plants and algae 
since it transforms electromagnetic energy into chemical energy. Due to its impor- 
tance in converting solar energy to a biologically utilizable form [1], the molecular 
mechanism of photosynthesis has been widely studied, especially in photosynthetic 
bacteria such as green sulphur bacteria, purple phototrophic bacteria, green filamen- 
tous bacteria, cyanobacteria, and heliobacteria [2-4]. 

The light harvesting process consists of several steps. First, the incoming photon 
excites the electronic state of pigments in the antenna complex of the photosystem. 
Then the excited state energy will be transported to the reaction center of the pho- 
tosystem. When enough excitation energy has been collected, a charge separation 
will be triggered by one or a pair of bactcriochlorophylls (BChls) at the reaction cen- 
ter; this will further induce an electron transfer via an electron transport chain and 
generate a proton gradient across the membrane. This proton gradient drives ATP 
synthase, which produces the most common energy form (ATP) utilized in organisms. 
Several reviews of the light harvesting process in purple bacteria are provided here 
[5-9]. 

In purple bacteria - a model organism in photosynthesis studies - the photosyn- 
thetic units are located on the inner membrane and are comprised of two types of 
light harvesting membrane pigment-protein complexes. One is light harvesting com- 
plex 2 (LH2), which is the peripheral antenna that mainly absorbs light. The other 
is light harvesting complex 1 (LHl), which surrounds a reaction center and forms the 
LHl-RC complex which acts as the core antenna. The main pigments that absorb 
light and contribute to the excitation energy transfer in these antenna complexes are 
carotenoid (Car) and bacteriochlorophyll (BChl). 

One of the most intriguing features of photosynthesis is the near-unity efficiency 
of energy conversion from light to charge [4] . The detailed mechanism of this delicate 
biological process remains elusive. The most well-studied systems are the Fenna- 
Matthews-Olson (FMO) complex of green sulfur bacteria and the LH2 complex of 
purple bacteria. The FMO complex transfers the exciton energy from the light har- 
vesting complex to the reaction center (RC), thus plays the role of an intermediate 
candidate which passes the exciton energy with an efficiency of almost 100% [10]. 
A long-lived quantum coherence was experimentally found in the pigment excited 
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electronic states [11], implying that quantum mechanical effects may be involved in 
the excitation energy transfer. Recently, Aspuru-Guzik et al. studied the effects on 
efficiency of the combination of quantum coherence and environment interaction via 
a non-Markovian approach developed from the Lindblad form [12-15]. Meanwhile, 
Ishizaki and coworkers [16, 17] reproduced the experimentally detected population 
oscillation time observed by Engel and coworkers [11] via the hierarchical equation 
of motion (HEOM) approach. The excitation energy transfer pathway in the FMO 
complex has been investigated by Skochdopole et al [18]. The role of pairwise entan- 
glement in FMO has also been investigated by Sarovar et al [19]. Recently we applied 
the scaled HEOM method developed by Shi et al. [20] to the study of the quantum 
evolution of the FMO complex, showing that the scaled HEOM method is compu- 
tationally more stable and accurate [21]. In addition, we also investigated the role 
of multipartite entanglement by direct computation of the convex roof optimization 
[22]. In addition, there are numerous theoretical studies on modelling the dynamics 
of excitation energy transfer in the FMO system [23-38]. 

For LH2 in purple bacteria, a long-time quantum coherence in the excited states 
of bacteriocholophylls (BChl) has also been found experimentally [5, 39-41], and the 
excitation energy transfer dynamics have been investigated theoretically [42, 43]. In 
this article we employed the scaled HEOM approach to study the quantum evolu- 
tion in the LH2 complex, and a quantitative description of the long-time quantum 
coherence is provided. The paper is organized as follows. In the method section, the 
theoretical framework of scaled HEOM as well as the coherence will be introduced. 
Next, the simulation result and the physical explanation for both population and 
coherence evolution will be investigated. 
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II. METHOD 



The system studied here is the BChls in the LH2 B850 ring of purple bacteria 
Rhodopseudomonas (Rps.) acidophila [44] (PDBID: INKZ). We employed the open 
quantum system method to explore the electronic excitation dynamics of this complex. 
The 18 B850 BChls of the LH2 ring are treated as the system, which couples to 
the environment modelled by an infinite set of harmonic oscillators. They interact 
with each other via system-environment coupling, and the total Hamiltonian can be 
described as: 

'H = T-is + "Hs + T-isBi (1) 

where l-Lsi and I-Lsb represent the Hamiltonian of the system, environment and 
system-environment coupling respectively. 

A. The Structure of the LH2 Complex 

The molecular structure of LH2 is highly symmetric, and usually possesses Cs 
or Cg symmetry depending on the different species of bacteria. For example, in 
Rhodospirillum {Rs.) molischianum [45] and Rps. acidophila [44, 46] LH2 has an 
eight-fold and nine-fold symmetry respectively. In Rps. acidophila each subunit is 
constructed by inner a and outer /3-peptides and containing 3 BChls and 1 Car. In 
the LH2 absorption spectrum, there are two significant peaks at 800 and 850 nm due 
to the Qy excitation of BChls. For each subunit there is one BChl that belongs to 
the 800 nm absorption peak and two BChls that contribute to the 850 nm peak. Due 
to symmetry, both of them form ring structures which are called the B800 and B850 
ring, respectively. 

In this paper, we will only study the quantum evolution of the B850 ring. We 
plan to explore the role of both rings (B800 and B850) during the energy transfer 
in the future. B850 ring is constituted by 18 BChls in nine transmembrane a, (3- 
polypeptide heterodimers and has a Cg symmetry [47]. The arrangement of 18 B850 
BChls is shown in Fig. 1. 

All 18 B850 BChls are marked clockwise, and the distance between BChls is de- 
fined by the distance between the Magnesium atoms corresponding to each BChl. 
The distance between adjacent B850 BChls in the same and different dimers are 
9.4 A and 9.1 A respectively [44]. The most popular theoretical models in this area 
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can be organized into three distinct groups, the Redfield equation, Forster theory 
and hierarchical equation of motion (HEOM). The Redfield equation [48, 49] assumes 
that the coupling between the system and environment is weak compared with the 
coupling within the system. Mathematically, we can treat the environment effects 
as a perturbation upon the system. The assumptions for the Forster theory are in 
total opposition to that of the Redfield equation [50, 51]. It assumes strong coupling 
between the system and environment so the interactions within the system can be 
neglected. As a result, the mathematical calculation can be simplified by treating 
the system Hamiltonian as a diagonal matrix. The HEOM does not require both as- 
sumptions. It can be applied to both strong and weak system-environment coupling, 
making it an accurate and comprehensive model [20, 52-54]. However, the numerical 
implementation requires a lot of computational resources. The electrostatic interac- 
tion between neighbouring B850 BChls is relatively strong due to the short distance 
between them, which makes the BChl electronic coupling at the same order as the 
electron-environment couphng. As a result, the HEOM is an ideal framework for this 
complex system. Our group has successfully employed the scaled HEOM approach 
to study the population beating and multiple particles entanglement evolution in the 
FMO complex [21, 22]. The simulation results also indicate that the scaled HEOM is 
helpful in reducing the required computational resources without losing any resolution 
or physical meaning. 

B. Hamiltonian of the system 

Under a single exciton basis set, the system Hamiltonian can be written as [55] 

TV 

ns = Y. \3) 01 + E d^') (^1 + 1^) 01) ' (2) 

where Sj represents the Qy excitation energy of the jth BChl and Jjk denotes the 
excitonic coupling between BChl j and k. There are only two different excitation 
energy for these 18 BChls. One is 12457,66 cm~^ and the other is 12 653,66 cm~^ 
[47]. They arrange alternatively 9 times [56]. In our simulation, we treated the energy 
of BChl 1 (SI) as 12 457,66 cm-^ and thus BChl 2 (S2) equals 12 653,66 cm-^ 

The non-nearest neighbor coupling between j and k was treated by dipole-dipole 
interaction, expressed as 
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Figure 1. (Color) The top view (from cytoplasmic to periplasmic side) of 2.0 A resolu- 
tion LH2 structure of Rhodopseudomonas acidophila. The 9 B800 and 18 B850 BChls are 
shown in green and cyan respectively, and the purple spheres represent Magnesium of the 
BChls. The 18 BChls are labelled clockwise from SI to S18, and the transition dipole of Qy 
excitation for each BChl is represented by red arrows. 



J _ \P\ Pj ■Pk-'i iPj ■ fjk) iPk ■ fjk) , . 



\rjk\ 



where pj and represents the unit vector of the j-th and k-th BChl dipole moment. 
\rjk\ indicates the distance between the Magnesium of the two BChls and fjk is the 
unit vector along fj^. The same dipole moment has been used for all BChls with the 
value of |pp/47reo set to 348 000 A^cm"^ [57]. For the two adjacent BChls in the same 
(SI and S2) and different dimers (SI and S18), the couplings are set as 363 cm~^ and 
320 cm~^ respectively [57, 58]. 
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C. The System-Environment Coupling 



The environment is described as a phonon bath, modelled by an infinite set of 
harmonic oscillators: 

AT AT NjB p2 

j=i j=i ^=1 ^« 

where rrij^, Uj^, Pj^, Xj^ are mass, frequency, momentum and the position operator of 
the ,^-th harmonic bath associate with the j-th BChl respectively. 

The Hamiltonian of the environment (T-Lb) and system-environment coupling 
{'Hsb) can be written as: 

N N 

'^SB = Y.Y. \^) 01 ^^-^ = E ^^^^^ (5) 

where Vj — \j) (j| and Fj = '^^Cj^Xj^. Cj^ represents the system-bath coupling 
constant between the j-th BChl and ^-th phonon mode. Here we assume each BChl 
is coupled to the environment independently. 

The initial condition can be written as: 

X(0)=p(0)®p^^ (6) 

where x, P, and are the total, system, and the bath density matrices, respectively. 
The system is in its ground state |0) (0| before being excited, and the environment 
is in a thermal equilibrium state. Once the system was excited by a laser pulse, 
the system suddenly changed to the specified state p(0), while the environment still 
resides at equilibrium. After t = 0, the system starts to evolve under the effect of the 
system-environment bath [59, 60]. The time evolution of the system can be obtained 
by tracing the environment p (t) = Tvb [x (t)] = Tr^ [e~'"*/'^x (0) e'^*/"] . 
The correlation function under phonon bath can be written as: 

oo 

C, {t) = -j <k.D, (a;) , (7) 

— oo 
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where ^ — I/ZcbT. We replaced the original spectral density (Eq. 8) with the Drude 
spectral density 



2A7 UJ 



(9) 



where A is the reorganization energy and 7 is the Drude decay constant. The corre- 
lation function can be written as: 



Cj{t > 0) = J2ck-e- 



(10) 



k=0 



with vq — 7, which is the Drude decay constant, Vk — 2k'!r/ f3h when k ^ 1 and Vk is 
known as the Matsuraba frequency. The constants Ck are given by: 



777 



cot 

2777 Vk 



for/c > 1 . 



(11) 
(12) 



Using the scaled HEOM approach developed by Shi et al [20] as well as the Ishizaki- 
Tanimura truncation scheme [52, 54], the time evolution of the density operator be- 
comes [21]: 



k 

N / 00 \ N K . — _ 

-E E 1^ ■K.K,p„l]-.i:Eyg(c.-v.P„--4P„-v,). (13) 

j=l \m=K+l •'"^Z 3=1 fe=0 V I "I 

where n is defined as one set of nonnegative integers n = {ni,n2,-- - ^un} — 
{{riio, rill • • • , riift-}, • • • , {nmP'Ni ■ ■ ■ , nNx}}- K is the truncation level of the cor- 
relation function and refers to change the value of rijk to rijk ± 1 in the global 
index n. The sum of Ujk is defined as tier Mc, Mc — Ylij k ''^jk- This is another trunca- 
tion in our numerical simulation. The density operator with all indices equal to is 
the system's reduced density operator (RDO) while all other operators are auxihary 
density operators (ADOs). 

Apart from the population evolution, the coherence evolution is also evaluated 
here. The coherence between site i and j is defined as \pij\ 
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D. Simulation Details 



During the numerical simulation, we tried different truncation levels for both the 
correlation function (K) and tiers (Afc) and checked the relative convergence. The 
highest truncation levels we can reach are K — 2 and N'c — 3. 

The reorganization energy (A) and Drude decay constant (7) are set at 200 cm~^ 
and 100 fs"^ respectively [61, 62]. The simulation at these truncation levels require 
at least 19.6 GB of memory, and for each integration step (1 fs) it takes about 70 s 
on a 2.3 GHz cpu core. 

Compared with the FMO complex [21], these truncation levels are relatively low 
as the system of the LH2 B850 BChls ring is much larger (18 x 18) than that of FMO 
situation (7 x 7). The time evolution of the LH2 system (Eq.(13)) was calculated by 
the 4th order Runge-Kutta method with a fixed time step of 1 fs. We calculated the 
excitation dynamics for up to 1.5 ps. 
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III. RESULTS AND DISCUSSION 



A. Exciton Population Dynamics 

In the non-dissipative case (Fig. 2(a)), the population at Si decreases rapidly and 
reaches its minimum value at 19 fs, while the populations of S2 and S18 continually 
increase and arrive at their maximum value of 0.334, 0.243 at 14, 13 fs respectively. 
When t — 36 fs, the populations of S5 and S15 become the maximally populated 
BChls. At 68 fs, the population of SIO dominates the whole LH2 system. Interest- 
ingly, the population of SI hits another maximum value of 0.410 at 128 fs, which 
is almost twice of the time for the population passed from SI to SIO (68 fs). This 
suggests that the excitation wave packet propagates both clockwise (towards S2) and 
counterclockwise (towards S18) and those two wave packets superpose at SIO and 
then continue traversing the entire ring. 

For the dissipative case (Fig. 2(b)), the populations change in a similar fashion to 
the non-dissipative situation for the first 50 fs. The simulation was run at tempera- 
ture T — 300 K. There is also population beating with the effect of the environment. 
The beating time of the population oscillation lasts roughly 150 fs. The populations 
of S2, S5 and SIO evolves to their first maximum value at 13, 35 and 66 fs respec- 
tively. Compared with the previous, non-dissipative case, these maxima are achieved 
simultaneously, which indicates the population oscillation of LII2 is due to the system 
itself rather than the environment. The environment plays a role in destroying the 
population beating, which is analogous to that in the FMO complex [21]. However, 
the population amplitude of each site and beating time are both smaller when com- 
pared to the FMO complex. This is reasonable since LH2 is much larger than the 
FMO complex. 

In contrast with the non-dissipative situation, the populations of all sites evolve 
to their equilibrium state during the time evolution. Furthermore, the difference 
for the peak arrival time also becomes larger as the evolution continues. It shows 
that the environment not only affects the population distribution but also facilitates 
the propagation speed of these exciton waves. The system arrives at the thermal 
equilibrium at 400 fs when T — 300 K. Due to the symmetry of the system, all odd- 
number sites will finally equilibrate to the same population and similarly with even 
sites. When two wave packets first reach SIO, the local maximal population is about 
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1.6 times that of the first local maximum of the S5 and S15 populations, and the time 
variance of this peak increases as well. Although the system is highly symmetric, the 
coupling between sites is non-trivial and makes the oscillation pattern complex, which 
is quite different than in the FMO system. Using another initial condition (S2 fully 
excited), we found that in the non-dissipative case (Fig. 3(a)), the time when two 
opposite wave packets superimpose is identical to the situation of SI initially excited. 

For the dissipative case (Fig. 3(b)), the amplitude of the first population maximum 
of adjacent sites corresponds to the strength of coupling. Since the coupling between 
BChls in the same dimer is stronger than that between 2 BChls of different dimers, 
larger population can be found in the BChl in the same dimer. When SI is initially 
excited, the population of SI decreases rapidly to about 0.1 and oscillates around 
this value for about 150 fs. However, when taking S2 as the initial state, it drops 
quickly to 0.2 at first and then decreases toward its equilibrium by about 200 fs. The 
population decrease of S2 is significantly larger than all other sites initially. However, 
the equilibrium population is found to be independent of initial conditions. 




Time / fs Time / fs 



Figure 2. (Color) The exciton population dynamics of B850 bacteriochlorophylls (BChls) 
with site 1 (SI) initially excited, (a) The population evolution of SI, S2, S5, SIO, S15, and 
S18 without dissipation (the system is isolated and uncoupled to bath), (b) The population 
dynamics of the same sites while the system is coupled to bath at room temperature T — 
300 K. The coherent energy transfer lasts about 150 fs and the whole system is equilibrated 
after 400 fs. The inset is a magnification of the first 250 fs dynamics. 
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Figure 3. (Color) The exciton population dynamics of B850 bacteriocliloropliylls (BClils) 
with site 2 (S2) initially excited, (a) The population evolution of SI, S2, S3, S6, Sll, and 
S16 without dissipation (the system is isolated and uncoupled to bath), (b) The population 
dynamics of the same sites while the system is coupled to bath at room temperature T — 
300 K. The coherent energy transfer lasts about 150 fs and the whole system is equilibrated 
after 400 fs. The inset is a magnification of the first 250 fs-dynamics. 



B. Coherence Dynamics 



From Fig. 4(a), it can be seen that the coherence lasted all the time for the non- 
dissipative case. The simulation result is valid as there are no exterior effects elimi- 
nating the coherence for an isolated system. Several important pairs are plotted in 
Fig. 4, which are S1-S2, S1-S5, Sl-SlO, S1-S18, S9-S10, and S9-S11. All the other 
pairs can be predicted through system symmetry. Although SI and SIO are well 
separated in distance (about 52.75 A), their coherence is still significant. 

On the other hand, when a heat bath (T = 300K) is coupled to the system (Fig. 
4(b)), the coherence drops quickly during the evolution. For a quantitative view of 
decoherence, the second peak of S1-S5 coherence at 300 K and non-dissipative case 
has been compared, and we found that it suppresses 39.27% at t = 27 fs when the 
system is dissipative. In addition, in the first peak of S9-S10 coherence (around 
t = 64 fs) the amplitude decreases 73.56% and clearly demonstrates the decoherence 
effect. At the nearly-equilibrium state, the significant coherence (close to 0.1) occurs 
between adjacent sites (S1-S2, S1-S18, and S9-S10) and between second nearest sites 
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(S9-S11). For the situation of S2 initially excited, the first peak of S2-S6 and the 1st 
peak of S9-S10 coherence show a decrease of 58.77% (at t = 26 fs) and 48.63% (at 
t = 56 fs) respectively in the dissipative case (Fig. 5(b)). This indicates that the 
environment plays an important role in reducing the coherence. 




50 100 150 200 250 50 100 150 200 250 



Time / fs Time / fs 

Figure 4. (Color) The first 250 fs-coherence dynamics of B850 bacteriochlorophyll under 300 
K with SI initially excited (a) with or (b) without coupling to a phonon bath. The coherence 
decays as the distance between the sites increases. There is a significant long-time coherence 
between site 1 and site 2 excitons even the system is equilibrated. The inset shows the 1 
ps-dynamics at the same condition. 

To obtain a holistic picture of the information contained in the density matrix, 
the absolute value of all density matrix elements have been plotted (Fig. 6). The 
diagonal elements represent the exciton population, and the off-diagonal elements are 
the coherence between two sites. The initial state of Fig. 6 is SI fully excited. After 
t = 0, the exciton wave packets propagate both clockwise (towards S2) and counter- 
clockwise (towards S18) and merge at SIO at t = 66 fs. Then, they traverse through 
the whole ring and back to SI (at ~124 fs). After that the system regains some coher- 
ence between adjacent sites. Even when the system is near its equilibrium (1.5 ps), 
there exists a significant coherence between several nearest neighbours (4 and 2 ad- 
jacent sites for odd and even number sites respectively) and thus provides additional 
evidence for the existence of long-lived coherence between adjacent chromophores in 
LH2 system. 
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Figure 5. (Color) The coherence dynamics of B850 bacteriochlorophyll under 300 K with 
S2 initially excited (a) with or (b) without coupling to a phonon bath. The inset shows the 
1 ps-dynamics at the same condition. 

C. Temperature Effects on Exciton Dynamics 

In order to investigate the effects of temperature on exciton dynamics, the popula- 
tion and coherence dynamics have been studied at cryogenic temperature (77 K), this 
being the temperature at which most 2D electronic spectroscopy experiments have 
been performed. 

The exciton population dynamics at 77 K are shown in Fig. 7. Similar to the 
300 K case, the exciton wave packets propagate both clockwise and counterclockwise. 
The time taken for two wave packets to merge at SIO is 67 fs, which is shorter than 
that of the isolated situation (68 fs), but longer compared with that for T = 300 K 
(66 fs). The time for the excitation wave packet returning to SI is 127 fs, which is 
also located between that for isolated and T = 3007^ cases. The same phenomenon 
can be observed by checking the time when S5 and SIO have their 3rd maximum 
(non-dissipative: 106 fs; 77 K: 105 fs; 300 K: 104 fs). The environment increases the 
population propagation speed and higher temperature further improves the speed. 

The population oscillation lasts over 300 fs at cryogenic temperature, which is 
almost twice that of at room temperature (~ 150 fs). More quantitatively, one could 
also check the 1st maximum population ratio of S5 (or S15) to SIO, which indicates 
the dissipative effect while the exciton wave packet travel from S5 to SIO. At non- 
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Figure 6. (Color) The time evolution of density matrices (absolute value) of the 18-sites 
system at 300 K. The initial condition is SI fully excited, (a) i = 35 fs, when S5 and S15 
have their maximum population, (b) i = 66 fs, when SIO hits its maximum population, (c) 
t = 124 fs, when both exciton wave packets traverse through the ring, (d) t = 1500 fs, the 
system is in the equilibrium state. 



dissipative, 77 K and 300 K, these ratios are 1.634, 1.106, and 0.651 respectively. 
Since the first maximum population appears at almost the same time for all three 
conditions, it is valid to say that the dissipation rate gets larger with increasing 
temperature. A consequence of this is that the system will take longer to equilibrate 
(~600 fs at 77 K) at low temperature. In other words, the system-environment 
coupling becomes stronger when the temperature increases; high temperature plays 
a positive role in driving the system into thermal equilibrium. 

The coherence dynamics at 77 K is shown in Fig. 8. The decay rate shows that the 
system decoheres slower compared to the case at 300 K. For the first 150 fs, all sites 
have relatively large coherence. On there other hand, the coherent energy transfer 
lasts about 300 fs, which is the same length as the population oscillation time. Again 
the second maximum of S1-S5 coherence (SI initially excited) and the first maximum 
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Figure 7. (Color) The excitation dynamics of B850 bacteriochlorophyll with (a) SI or (b) S2 
initially excited under 77 K. The coherent transfer between exciton lasts about 300 fs and 
the whole system is equilibrated after 600 fs. Longer population oscillation time is observed 
comparing to 300 K. Also around 127 fs (SI initially excited) and 125 fs (S2 initially excited) 
the initial excited site has a local maximum with a larger width, which indicates when the 
two opposite exciton packets superpose. The inset is a magnification of the first 250 fs- 
dynamics 



of S2-S6 coherence (S2 initially excited) have been compared to the non-dissipative 
case, and in both cases the coherence only be reduced by 18.64% at t = 28 fs and 
32.87% at t = 29 fs, respectively. These results show that the dissipative effect occurs 
slower than the high temperature case. We conclude that the temperature plays a 
negative role in coherence conservation. Lower temperature is helpful for the system 
to maintain the coherence. 

After checking the absolute value of density matrix elements at 77 K (Fig. 9) 
starting from SI excitation condition, we found that when S5 and S15 have their 
maximum, the S5-S15 coherence is significantly larger than that at room temperature. 
The same also occurs when SIO hits its first maximum, significant coherence appears 
between SIO and odd number sites all over the ring. This phenomenon is more 
apparent at low temperature. When the two opposite exciton wave packets propagate 
back to SI, in spite of a larger SI population at low temperature, the coherence 
between S1-S7 and S10-S13 is also significant. For the final equilibrium state (at 1.5 
ps) at 77 K, while the even number sites have coherence mainly with the nearest two 
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Figure 8. (Color) The coherence dynamics of B850 bacteriochlorophyll under 77 K with (a) 
SI or (b) S2 initially excited. The oscillations of the coherence last about 250 fs, which 
is significantly longer than in the 300 K. The coherence between sites also decays as the 
distance increases. The inset shows the 1 ps-dynamics at the same condition. 

neighbours, the odd number sites have coherence with its four nearest neighbours and 
the other odd number sites. This effect is absent in the high temperature case. 



IV. CONCLUSIONS 

In summary, we employed the scaled hierarchical equation of motion (HEOM), 
which has been proved to be computationally more efficient [21], in studying the pop- 
ulation and the coherence evolution of the LH2 B850 for both cryogenic (T = 77 K) 
and physiological (T = 300 K) temperatures. Moreover, the short-time dynamics 
and the thermal equilibration of excitation energy transfer within the ring are sys- 
tematically investigated by visualizing the reduced density operator. 

Similar to the situation for the FMO complex, there are both population and co- 
herence oscillations during the time evolution in the LH2 system. In addition, the 
oscillation time for population and coherence is almost the same for the same system. 
High temperatures shorten the coherence time and suppress both oscillation ampli- 
tude of coherence and population. However, the beating time for LH2 can last 150 
fs and 300 fs at physiological and cryogenic temperatures, respectively. This is much 
shorter compared with the FMO complex (coherence oscillation can last 400 and 650 
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Figure 9. (Color) The time evolution of density matrices (absolute value) of the 18-sites 
system at 77 K. The state with fully excited SI is used as the initial condition, (a) t = 36 
fs, when S5 and S15 have their maximum population, (b) t = 67 fs, when SIO has its 
maximum population, (c) t = 127 fs, when both exciton wave packets traverse through 
about the whole ring, (d) t = 1500 fs, the system is in a nearly equilibrated state. 



fs at the same temperature). Besides this, the exciton packet transfer both clockwise 
and counterclockwise during the evolution for LH2 system. The period for exciton 
transferred is 128, 127 and 124 fs for the non-dissipative, cryogenic and physiolog- 
ical temperature situations; this indicates that the high temperature increases the 
propagation speed of the exciton wave packet. Moreover, although the coherence os- 
cillation vanishes rapidly, the coherence between adjacent BChls revives and remains 
in the system for a long time, which is consistent with the result of 2D electronic 
spectroscopy from other experimental studies. This phenomenon may give us some 
clues about coherence could be maintained in such a noisy environment and may 
shed a light on designing a quantum system which can preserve coherence at room 
temperature. 

Recently, the optimized HEOM coupled with the Drude spectral density has been 



18 



proposed [63, 64], which improves the efficiency in numerical calculations. It will be 
our future task to apply this approach to the current model for simulating larger 
systems, such as the LH2 B800-B850. 
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